Asymmetric vortex solitons in nonlinear periodic lattices 
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We reveal the existence of asymmetric vortex solitons in ideally symmetric periodic lattices, and 
show how such nonlinear localized structures describing elementary circular flows can be analyzed 
systematically using the energy-balance relations. We present the examples of rhomboid, rectangular, 
and triangular vortex solitons on a square lattice, and also describe novel coherent states where the 
populations of clockwise and ant i- clockwise vortex modes change periodically due to a nonlinearity- 
induced momentum exchange through the lattice. Asymmetric vortex solitons are expected to 
exist in different nonlinear lattice systems including optically-induced photonic lattices, nonlinear 
photonic crystals, and Bose-Einstein condensates in optical lattices. 
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Vortices are fundamental objects which appear in 
many branches of physics [l|. I n optics, vortices are as- 
sociated with phase dislocations (or phase singularities) 
carried by diffracting optical beams |2|, and they share 
many common properties with the vortices observed in 
superfluids and Bose-Einstein condensates (BEC) Q. 
Vortices may emerge spontaneously or can be generated 
by various experimental techniques, and they are impor- 
tant objects for fundamental studies also having many 
practical applications. Nonlinear self-action effects are 
inherent to many physical systems, and they support the 
existence of localized vortex-like structures in the form of 
vortex solitons when the flow around the singularity does 
not change the density distribution, the property natu- 
rally observed for axially symmetric structures in homo- 
geneous media [Fig.Qfa)]. 

Periodic lattices, such as photonic structures for laser 
beams or optical lattices for atomic BECs, allow for a 
strong modification of the wave propagation, which also 
depends on the energy density. Recently, it was demon- 
strated experimentally that photonic lattices can support 
stable off-site (small radius) and on-site (larger radius) 
vortices on a square lattice |4| . The square-like profiles of 
such vortices resemble a homogeneous vortex modulated 
by the underlying periodic structure [Fig. Qfb)]. Vortex 
solitons describe elementary circular energy flows, and 
therefore their properties are intrinsically linked to the 
wave transport mechanisms in the underlying periodic 
lattice. In this Letter, we demonstrate that this con- 
nection is highly nontrivial due to an interplay between 
nonlinearity and periodicity. We predict analytically and 
confirm by numerical simulations that even ideally sym- 
metric periodic structures can support robust asymmet- 
ric vortex solitons, and our approach allows for a sys- 
tematic study of such novel types of singular states. We 
find that some symmetries are always allowed, whereas 
other configurations, such as a triangular vortex shown 
in Fig. [He), may exist under certain conditions derived 
from a balance of the energy flows. Additionally, we pre- 
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FIG. 1: Schematic of (a) conventional vortex in a homoge- 
neous medium, (b,c) symmetric and asymmetric vortex soli- 
tons in a periodic lattice. Arrows indicate the direction of the 
phase change and associated energy flow around the screw 
dislocation (phase singularity) marked at the vortex center. 



diet the existence of fully coherent states in the form of 
vortices which exhibit a charge flipping, i.e. a periodic re- 
versal of the energy flow. This effect can occur for larger- 
radius vortices, uncovering a key difference between the 
on-site and off-site vortex states. 

We consider the nonlinear propagation of an optical 
beam in a two-dimensional periodic lattice 5] described 
by the dimensionless nonlinear equation, 
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y + V(x,y)y-g(x,y,\y\ 2 )^ = 0, (1) 



where stands for the transverse Laplacian, ^(x, y, z) 
is the complex field amplitude, x, y are the transverse 
coordinates, z is the propagation coordinate, D is the 
diffraction (or dispersion) coefficient. Function V defines 
a periodic potential of the two-dimensional lattice, and 
the function Q characterizes a nonlinear response. Sim- 
ilar mathematical models appear for describing the self- 
action effects in nonlinear photonic crystals 6], and the 
nonlinear dynamics of atomic BEC in optical lattices • 
Nonlinearity in Eq. can compensate for the 

diffraction-induced beam spreading in the transverse di- 
mensions leading to the formation of stationary struc- 
tures in the form of spatial solitons, ty(x,y,z) = 
ip(x, y)e l $ z , where x/>(x,y) is the soliton envelope, and 
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P is a nonlinear shift of the propagation constant, the 
soliton parameter. Periodic lattices can modify dramat- 
ically the soliton properties, as was demonstrated for 
waveguide arrays and dynamically-induced periodic pho- 
tonic lattices 5j. In particular, periodic lattices can sta- 
bilize and suppor t propagation of discrete vortex soli- 
tons || El ll(t fi"H which have recently been observed in 
experiment |4|. 

In order to analyze the vortex-like structures in a pe- 
riodic potential, we write the field envelope in the form, 
i/j(x, y) = \i/j(x, y) \ exp[z(^(x, y)], and assume that the ac- 
cumulation of the phase ip around a singular point (at 
i/j = 0) is 27rM, where M is an integer topological charge. 
We consider spatially localized structures in the form of 
vortex-like bright solitons with the envelopes vanishing 
at infinity. Such structures can exist when the soliton 
parameter (3 is placed inside a gap of the linear Floquet- 
Bloch spectrum of the periodic structure 0, Q • 

The profiles of four-lobe symmetric vortex solitons dis- 
cussed so far in both the theory || 0, [HI E| and ex- 
periment Q resemble closely a ring-like structure of the 
soliton clusters m homogeneous media. Then, we 
look for novel vortex solitons of arbitrary symmetry as a 
superposition of a number of the fundamental solitons, 

N-1 

V(x, y-z)=J2 A n (z)i/j s (x -x n ,y- y n )e i(3z , (2) 

n=0 

where ip s are the profiles of the individual fundamental 
solitons, n = 0, . . . , N — 1, N is the total number of soli- 
tons, (x n ,y n ) are the soliton positions, and A n is the 
scaling coefficient defining the phase of the n-th soliton 
and variation of its amplitude due to interaction with 
other solitons. In contrast to the case of a homogeneous 
medium Q, the positions of the individual solitons are 
fixed by the lattice potential, provided the lattice poten- 
tial is sufficiently strong. In order to determine the soli- 
ton amplitudes, we present Eq. JU in the Hamiltonian 
form, idV/dz = SH(x,y,9 , **) , and derive, after 
substituting Eq. (J2J) into the full Hamiltonian, the re- 
duced Hamiltonian H 3 (x,y,A n ,A^) . The resulting am- 
plitude equations can be written in the form, 

^~dz~ ~ 5 A* ~ ~ ^ 1 Cnrn ^ rn ~ 0(\A n \ )A n — F ni (3) 

where c nrn = c mn = C nrn /C nn , C nm = j j ips(x 
x n ,y - y n )i/j*(x - Xm,y - y m )dxdy are the coupling coef- 
ficients, G is the effective local nonlinearity, and F n de- 
fines the nonlinear coupling terms such as ~ A mi A^ Tl2 A n . 
We note that the approximation (|2j) is valid when 
Cn+ rp., G, F /A n ~ e <C 1. As was demonstrated ear- 
lier PJJ [fl > under such conditions the amplitudes A n are 
only slightly perturbed, and we can seek stationary solu- 
tions of Eq. (|3J) corresponding to vortex solitons by means 
of the perturbation theory: A n = [1 + 0(e)] exp[i(p n + 



iO(e)). In Eq. (|3J), the nonlinear coupling terms are pro- 
portional to the forth-order overlap integrals and, there- 
fore, F n <C s 2l rn c nm A rn . Then, in the first order we 
obtain a general constraint on the soliton phases </? n , 

N-1 

^ Cnrn Sin((^ m - lf n ) = 0. (4) 
m=0 

In the sum Q , each term defines the energy flow between 
the solitons with numbers n and m, so that the equa- 
tions (HJ represent a condition for a balance of the energy 
flows which is required for stable propagation of a soliton 
cluster and the vortex-soliton formation. These condi- 
tions are satisfied trivially when all the solitons are in- or 
out-of-phase. We note that Eq. ((SJ with F n = has the 
form of a discrete self-trayyinq equation, which appears 
in different physical contexts However, nontriv- 

ial solutions of Eqs. (@J) corresponding to the vortex-like 
soliton clusters have been analyzed only for symmetric 
configurations, and even then some important solutions 
have been overlooked, as we demonstrate below. More- 
over, we show that the existence domains of asymmetric 
vortex-like solutions are highly nontrivial, due to specific 
properties of the coupling coefficients calculated for real- 
istic periodic lattices. 

In order to provide a direct link to the recent ex- 
perimental results 4], first we apply our general ana- 
lytical approach to describe the vortex solitons in two- 
dimensional optically-induced lattices created in a pho- 
torefractive crystal. Then, the diffraction coefficient in 
Eq. is defined as D = z^X/ (Aimsx^), where xo and 
zo are the characteristic length-scales in the transverse 
and longitudinal spatial dimensions, respectively, no is 
the average medium refractive index, and A is the vac- 
uum wavelength. The lattice potential and nonlinear 
beam self-action effect are both due to the photorefrac- 
tive screening nonlinearity, 

g-v = 7 {i b + i 0S ^(^) s ^(^-) + \^y\ 

(5) 

where 7 is proportional to the external bias field, 1^ is 
the dark irradiance, and io is the intensity of interfering 
beams that induce a square lattice with the period d (see 
details in Refs. HQ). 

We use nonlinear localized (soliton) solutions obtained 
by using the ansatz (|2j) as initial conditions for a numer- 
ical relaxation algorithm applied to the model Q,©- 
We start by considering a four-soliton vortex on a square 
lattice [Fig. Efa)]. The previously reported solutions 
[Figs. E[b,c)] have both the reflection and 90° rotational 
symmetries, similar to the underlying lattice. Accord- 
ingly, coi = C12 = c 2 3 = c 30 and c 02 = C13, and Eqs. (@J) 
have a nontrivial-phase solution with cp n = nir/2 corre- 
sponding to a charge- one symmetric vortex. However, 
our analytical approach allows us to predict more gen- 
eral, reduced- symmetry vortex solitons when the 90° ro- 
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FIG. 2: Examples of the vortex- type soliton structures with 
various symmetries in a square lattice potential shown in (a) . 
(b,c) Off-site and on-site symmetric square vortex solitons 
(7 = 8,/? = 5.5); (d) rhomboid configuration (7 = 5, (3 = 3.5) 
with a topological charge +1. Shown are the intensity profiles 
(top) and phase structures (bottom). 

tational symmetry is removed. In particular, we find a 
rhomboid configuration that has a vortex charge with 
cp n = nir/2 [Fig. Efd)]. We found that rhomboid vortices 
are remarkably robust suggesting that their generation in 
experiment can be possible by using elliptically shaped 
singular beams. 

Most remarkably, for both square and rhomboid con- 
figurations, we find that the balance equations (@J admit 
more general exact solutions, which were overlooked in 
the earlier studies |l5| of Eq. ((SJ, namely cpo = if 2 — tt 
and if 1 — cps — 7r, where the phase difference cpi — cpo 
is arbitrary. These novel solutions describe a family of 
vortex solitons having the same intensity profile but dif- 
ferent phase structure. Due to such a degeneracy, a small 
change in the amplitude of two opposing solitons can ini- 
tiate a slow variation of the free phase, f\ — cpo ~ nz\ 
this regime corresponds to a periodic flipping of the vor- 
tex charge. Although the general equations are only 
satisfied when ipi — cpo = 0, 7r/2, . . ., the charge-flipping ef- 
fect can be induced by a finite perturbation when the the 
nonlinear coupling terms (F n ) are small. Indeed, we find 
that a closely packed vortex-like state shown in Fig.^b) 
is resistant to the charge flipping effect, similar to vortices 
in single- well potentials whereas the vortex shown 
in Fig. |2f c) can exhibit the charge flipping after increas- 
ing the amplitudes of two opposing solitons by 7%, as 
shown in Fig. 03 The bottom plot clearly illustrates that 
the solitons strongly exchange energy, however the flows 
always remain balanced. These are novel coherent states, 
where the populations of clockwise and anti-clockwise ro- 
tational modes change periodically due to nonlinearity- 
induced momentum exchange through the lattice. 

We now consider another example, vortices with rect- 
angular arrangement of the fundamental solitons. In this 
case, coi = C23, C03 = C12, and C02 = £13- With no loss 
of generality, we assign the soliton numbers such that 
coi < C03, and then nontrivial solutions of Eqs. (0} are: 

(Pl-(f0= COS -1 (-C 2/C03), f2~f0 = COS -1 (-Cqi / Cqs) , 
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FIG. 3: Charge flipping effect for a square vortex of Fig. Etc) 
induced by a 7% increase of the amplitudes of two opposite 
solitons. Top: Snapshots at increasing propagation distances 
showing the unchanged intensity profile. Arrows show the 
energy flow. Bottom: energy flows between the solitons char- 
acterized by sinusoidal functions of the phase differences be- 
tween the opposite (dashed line) and neighboring (solid and 
dotted lines) solitons according to Eq. Q. 




FIG. 4: Examples of asymmetric vortex solitons with rect- 
angle configurations. Top: intensity profiles with arrow indi- 
cating the energy flow; Bottom: phase structures. 



and cps + (fo = ipi + ^2- Since an inverse cosine function 
has two branches, there exist four solutions correspond- 
ing to two pairs of positively and negatively charged vor- 
tices with different positions of singularity. We find that 
the singularity is always shifted away from the center 
of rectangle along its longer dimension, as illustrated in 
Figs. H{a,b). This happens due to a highly asymmetric 
phase structure of the vortex, which in turn can lead to 
deformations of the vortex intensity profile resulting in a 
trapezoid-like shape shown in Fig.0fa). 
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FIG. 5: (a) Existence region in the plane of normalized 
coupling coefficients for the three-soliton vortex; (b) Exam- 
ple of a vortex soliton with the phase configuration ipo = 0, 
cpo ~ 0.9817T, Lp2 — — 0.50977T, supported by a self-focusing 
photorefractive screening nonlinearity (7 = 8, Iq = 1, D = 1, 
(3 = 5.5); (c) Dependence of the coupling coefficients on the 
soliton propagation constant (3 inside the band-gap for a con- 
figuration show in (d); (d) Example of a vortex soliton with 
the phase configuration ipo = 0, (pi ~ 0.6247T, <p2 — — 0.6257T, 
in repulsive atomic condensate (7 — 1, Vb — 2.5, D — 1) 
corresponding to marked points in plots (a,c) (/3 = 2.5). 

Finally, we mention another remarkable example of 
asymmetric vortex solitons in the form of a triangular 
structure. We find that Eqs. possess nontrivial so- 
lutions for N = 3 only when the normalized coupling 
coefficients, c\ = C12/C01 and C2 = C02/C01, satisfy the fol- 
lowing conditions, (1 + |ci| -1 ) -1 < |c2 1 < |1 — 
and we show the existence domain in Fig. EJa). It fol- 
lows that highly asymmetric triangular vortices are not 
allowed. On the other hand, isosceles configurations are 
always possible if the soliton base is sufficiently narrow, 
so that condition c\ = C2 < 2 is satisfied. We show an 
example of such vortex soliton in Fig.[Sfb). 

In order to underline the generality and broad applica- 
bility of our results, we also search for asymmetric vor- 
tex solitons in repulsive atomic Bose-Einstein conden- 
sates loaded onto a two-dimensional optical lattice. In 
the mean-field approximation, this system is described 
by the Gross-Pitaevskii equation (Q with 

V = Vo[sm(7rx/d) + sin(7ry/d)]; G = 7l^| 2 , (6) 

where 7 > 0, and Vb is the depth of an optical-lattice 
potential. Repulsive nonlinearity can lead to wave lo- 
calization in the form of gap solitons 7], which exhibit 
effective long-range interaction. This allows, in particu- 
lar, the existence of triangular vortices with a wider base 



[Fig. Etd)] compared to the conventional solitons in the 
self- focusing regime [cf. Fig.EJb)]. We find that such tri- 
angular gap vortices exist throughout the whole band-gap 
since c\ = C2 < 2 [Fig. Etc)], except for the immediate 
vicinity of the spectrum band edges where our approxi- 
mation is not applicable. 

In addition, we have verified that all vortex-like struc- 
tures discussed above are created by strongly coupled soli- 
tons. This is an important issue, since it may be formally 
possible to construct various composite states of almost 
decoupled solitons, which exchange the energy at an ex- 
tremely slow rate, so that Eqs. (HJ) are nearly satisfied. 
In contrast, in the cases discussed above there is a strong 
internal energy flow within the phase-locked and charge- 
flipping vortices, as seen in Fig. to Fig. El 

In conclusion, we have revealed that periodic lattices 
can support different types of robust asymmetric vortex- 
like nonlinear localized structures. Such vortices resem- 
ble the soliton clusters trapped by the lattice, but they 
are associated with a non-trivial power flow. We have 
presented the examples of novel vortex solitons on a 
square lattice, and other solutions can be obtained and 
analyzed using the general energy-balance relations. We 
believe our findings will initiate the experimental efforts 
to observe such vortices in optically-induced photonic 
structures, Bose-Einstein condensates in optical lattices, 
photonic crystals, and photonic crystal fibers. 
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